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We study non-linear adiabatic connection paths in density-functional theory using modified 
electron-electron interactions that perform a long-range/short-range separation of the Coulomb in- 
teraction. These adiabatic connections allows to define short-range exchange-correlation potentials 
and short-range local exchange-correlation energies per particle that we have calculated accurately 
for the He and Be atoms and compared to the corresponding quantities in the local density approx- 
imation (LDA). The results confirm that the LDA better describes exchange-correlation potentials 
and local exchange-correlation energies per particle when the range of the interaction is reduced. 
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I. INTRODUCTION 

The adiabatic connection procedure of density func- 
tional theory (DFT) Q connects the fictitious non- 
interacting Kohn-Sham (KS) Q system to the physical 
one by continuously switching on the electron-electron in- 
teraction. The connection is carried out at constant den- 
sity and is controlled by a interaction parameter \i. The 
partially interacting TV-electron system for a given value 
of [i along this connection is described by the Hamilto- 



H» = T + VJi + V^, 



(1) 



where T is the kinetic energy operator, V£, = 
is the modified electron-electron interac- 



<j "eeV y ) 

tion operator and V>* = J2i v>1 ( r i) ^ s the corresponding 
local external potential which ensures that this fictitious 
system has the same ground-state density n than that of 
the physical system. The ground-state wave function of 
this partially interacting system is denoted by \&' x . 

In the more common adiabatic connection (see, e.g., 
Ref.0), the electron-electron interaction is switched on 
linearly. However, other non-linear paths are possi- 
ble Q. In particular, in view of constructing a multi- 
determinantal extension of the KS scheme of DFT, it has 
been proposed 0, IS 0> IS 01 \Ei EJ to choose a modified 
interaction v£ e (r) representing the long-range part of the 
Coulomb interaction. In previous works 0, IS IS El , this 
long-range interaction has been defined by 



<e,crf( r ) 



erf (fir) 



(2) 



referred to as the erf interaction. More recently [l2| . 
a sharper long-range/short-range separation has been 
achieved by taking for the long-range interaction 
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erf (cj-ir) 
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(3) 
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referred to as the erfgau interaction. This modified in- 
teraction was first introduced by Gill and Adamson [l3T ] 

1/2 

in another context. In Eq. J^Jl, c = (l + 6^3) ~ 3.375 
is a constant chosen so that 1/fi roughly represents the 
range of the modified interaction as for the erf interac- 
tion [13. For both the erf and erfgau interactions, the 
modified interaction v£ e vanishes at fj, = and reduces 
to the full Coulomb interaction 1/r as fi — > 00. Conse- 
quently, the Hamiltonian of Eq. is the KS non- 
interacting Hamiltonian at fi — and the physical Hamil- 
tonian as /i — ► 00 . 

With these modified interactions, Eq. defines a fic- 
titious system with long-range interactions from which 
one can define a universal long-range density functional 



(4) 



and its short-range complement F M = F — F^ where F 
is the usual Coulombic universal functional. This short- 
range universal functional F 11 can be decomposed as 



pM — JJP 



(5) 



where = 1/2 JJ n(ri)n(ri)v£ e (ri2)dridr2 is the 

complement short-range Hartree energy functional and 
E£ c is the unknown complement short-range exchange- 
correlation energy functional. Finally, the ground-state 
energy of the physical system is given by 



v ne {v)n{v)dv + U^ + E^ c , (6) 



where v ne (v) is the nuclei-electron potential. Accord- 
ingly, the potential appearing in Eq. Q can be de- 
composed as 



^(r)=u ne (r)+<(r)+< c (r), 



(7) 



where v^(r) — SU^ /Sn(r) is the complement short-range 
Hartree potential and v% c (r) = 5Eg c /Sn(r) is the com- 
plement short-range exchange-correlation potential. As 
in the usual Kohn-Sham case, the short-range exchange- 
correlation energy can be decomposed into exchange and 
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correlation contributions, Eg c — Eg + E£. Correspond- 
ingly, the short-range exchange-correlation potential is 
decomposed as vg c (r) = vg(r) + vg(r). 

In this approach, the central quantity to approximate 
is the short-range exchange-correlation energy functional 
Eg c . Once an approximation has been chosen for Eg c , the 
potential appearing in Eq. is deduced by Eq. J7J), 
the multi-determinantal wave function 'E* is calculated 
and the ground-state energy E is deduced by Eq. JBJ. 
Previous applications of the method to small atomic and 
molecular systems 0, H, 0] show that, for a reasonable 
value of n, (^^{T + V^S 11 ) can be efficiently approxi- 
mated by standard wave function methods using a few- 
determinantal wave function 4'^. 

In practice, the short-range exchange-correlation en- 
ergy is conveniently expressed as 



W xc = \ n(v)e» c (v)dv. 



(8) 



where e^ c (r) is a short-range local exchange-correlation 
energy per particle that is not uniquely defined. The 
local density approximation (LDA) for e% c is 



? fi,LT>A 



(r) 



F /j,unif 



(n(r)), 



(9) 



where e^ un (n) is the short-range exchange-correlation 
energy per particle of a uniform electron gas with modi- 
fied interaction 

HEl. The LDA short -range exchange- 
correlation potential is 



<i LDA (r) 



'd( 



n e 



(n))' 



dn 



(10) 
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The performance of LDA for Eg c has previously be in- 
vestigated [H, |jj 0, Q . It has been found in particular 
on atomic systems that the LDA is very accurate for 
sufficiently large /i [T^|. For example, for the He atom 
with \x > 2, the LDA error on both Eg and E£ is less 
than 1 mil. However, good accuracy on integrated quan- 
tities does not always ensure a similar accuracy on local 
quantities. In this work, we analyze further the LDA for 
short-range exchange-correlation effects by performing a 
detailed local analysis. Indeed, for both the erf and er- 
fgau interactions, we have computed for the He and Be 
atoms accurate exchange-correlation potentials v% c and 
accurate local exchange-correlation energies per particle 
e£ c in two possible definitions and compared to the cor- 
responding LDA quantities. 

The paper is organized as follows. In Sec.[Hj the calcu- 
lation of accurate potentials v M is explained and the cor- 
responding exchange-correlation potentials are compared 
to the LDA ones. In Sec. IIIII accurate local exchange- 
correlation energies per particle are calculated in two 
possible definitions and compared to the LDA local en- 
ergy per particle. Finally, Sec. IIVI contains concluding 
remarks. 

Atomic units will be used throughout this work. 



II. SHORT-RANGE 
EXCHANGE-CORRELATION POTENTIALS 

The method we use for calculating accurate external 
potentials by optimization is explained in Refs.lldandllOl 
We rapidly recall here the principal aspects. 

The potentials i> M (r) are determined by using the 
Legendre transform formulation of the universal func- 
tional mm 



(11) 



F»[n] = sup^E^] - J ^(r)n(r)dr|, 



where E^lv^] is the ground-state energy of the Hamil- 

tonian T + V£, + X^ M ( r i)- If n is chosen to be the 
physical density of the system, and if n is assumed to 
be w-representable in the presence of the interaction V^ e , 
then the supremum in Eq. (|llfl is a maximum reached 
for the desired w M and . In practice, an accurate 
density n is computed by multi-reference configuration- 
interaction calculation with single and double excitations 
(MRCISD) [HIU and the potential to optimize v^(r) 
for atomic systems is expanded as 



n 



v"-(r) y^''- 
i=l 



c 

r 



(12) 



where Cj are the optimized coefficients, pi are some 
fixed integers (—1 or 2), 7, are fixed exponents cho- 
sen so as to form an even-tempered basis set (typically, 
7i G [10 -3 , 5.10 4 ]), and C is a constant enforces the cor- 
rect asymptotic behavior for r — > 00. For the Kohn-Sham 
case (/x = 0), this asymptotic behavior is determined by 
the nuclei-electron, Hartree and exchange contribution 
to the potential, giving C = —Z + N — 1 (N and Z are 
the electron number and nuclear charge, respectively). 
For finite /i, the short-range Hartree and exchange po- 
tentials are exponentially decreasing at infinity, so that 
it remains only the nuclei-electron contribution, giving 
C = —Z. Actually, the maximizing potential v M (r) in 
Eq. i|ll|) is defined only up to an additive constant. The 
potential of Eq. I|12|) is the one which goes to as r — > 00. 

The maximization of Eq. Qllfl is carried out with the 
Simplex method |2l|. For a given potential, E^[v^] 
is computed at MRCISD level using the Molpro pro- 
gram |22| with modified two-electron integrals (see Ap- 
pendix A of Ref . IT2I) . Beside the asymptotic behavior for 
r — > 00, v^(r) ~ C/r, the behavior of the potential at 
the nucleus r = 0, u M (r) ~ —Z/r, is also imposed dur- 
ing the optimization. Large one-electron even-tempered 
Gaussian basis sets are used for all systems (see Refs.ll6l 
and [l(J for more details) . The quality of the obtained 
potential is assessed using the Zhao-Parr criterion |2jj 

A= l_ /•/•(fi(r 1 )-n(r 1 ))(n(r 2 )-n(r a) ) 
2 J J |ri-r 2 | 

where h is the density given by the approximate potential 
. The largest values of A obtained with our potentials 
are of order 10~ 8 . 
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FIG. 1: Exchange potentials ri?£(r) for He with the erf inter- 
action with fj, = 0, 0.5 and 2: the potentials obtained in the 
present work (full curves) are compared to the LDA potentials 
[Eq. fWl. dashed curves . 



For two-electron systems, the exchange-correlation po- 
tential can be easily decomposed into exchange and cor- 
relation contributions. Indeed, the exchange potential is 
known from the Hartree potential, vtj:(r) = —v£{r)/2, 
and the correlation potential is simply obtained by dif- 
ference v£(r) — v£ c (r) — v£(r). Figure Q] shows the ex- 
change potentials of the He atom with the erf interaction 
for fi = 0, 0.5 and 2, together with the corresponding 
LDA potentials. To better visualize the asymptotic be- 
havior as r — > oo, rv(£(r) is plotted instead of v£(r). In 
the KS case (/i = 0), the LDA exchange potential largely 
differs from the accurate one. In particular, the LDA 
potential has not the correct Coulombic asymptotic be- 
havior u£ =0 (r) ~ — 1/r as r — ► oo. When n increases, 
the range of the accurate potentials decreases and rv£ (r) 
goes to as r — > oo. When \x increases, the asymptote is 
reached at smaller and smaller values of r. The important 
observation is that the accuracy of the LDA potentials 
increases with /x. At /i — 2, at the scale of the plot, the 
LDA potential is nearly identical to the accurate one. 

The accurate correlation potential for the He atom 
with the erf interaction for fi = 0, 1 and 3 is reported 
in figure[21 together with the LDA correlation potentials. 
In the KS case (/i = 0), the correlation potential calcu- 
lated by Umrigar and Gonze |24| by inversion of the KS 
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FIG. 2: Correlation potentials v£ (r) for He with the erf in- 
teraction with jj, = 0, 1 and 3: the potentials obtained in the 
present work (full curves) are compared to the LDA potentials 
[Eq. JTIJJ, dashed curves] . For fi = 0, the correlation potential 
calculated by Umrigar and Gonze |24| (dotted curve) is also 
shown. 



equation with a very accurate density is also reported. 
This potential agrees well with our potential except at 
very small r where the precision of our calculations (us- 
ing Gaussian basis sets) does not allow to extract the 
correlation potential. The LDA correlation potential for 
fx = is a poor approximation to the accurate, struc- 
tured potential. As for the exchange potential, the range 
of correlation potential is reduced when \x increases. At 
large the correlation potential has less structure which 
enables the LDA to perform better but on average only. 

The potentials obtained for the Be atom with the 
erf interaction are plotted in figure |3J At fi = 0, 
the exchange-correlation potential calculated by Umri- 
gar and Gonze (25J using quantum Monte Carlo methods 
is also reported. In comparison to the He atom, the accu- 
rate potentials for small /i now exhibit a shell structure 
with a jump separating the core region (r < 1) and the 
valence region (r > 1). The LDA potentials does not re- 
produce well this jump. When \x is large enough (/i > 1) 
so that the valence region is cut off, the shell structure 
disappears in the potential and the quality of the LDA 
potentials is improved. 

The potentials with the erfgau interaction are quali- 
tatively similar to the potentials with the erf interaction 
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FIG. 3: Exchange-correlation potentials rv£ c (r) for Be with 
the erf interaction with fi = 0, 0.5 and 1.5: the potentials 
obtained in the present work (full curves) are compared to 
the LDA potentials [Eq. mun . dashed curves]. For [x = 0, 
the exchange- correlation potential calculated by Umrigar and 
Gonze |25l| (dotted curve) is also shown. 



A. Short-range local exchange-correlation energies 
per particle from pair densities 

It is natural to define a short-range local exchange en- 
ergy per particle by 

^ ,P< Vi) = \J ^(ri,r 2 )<(r 12 )dr 2 . (14) 

where n x (r\,T2) is the exchange hole calculated from 
the one-particle KS density matrix: n x (ri,r2) = 
— |n5 <s (ri, r2)| 2 /(2n(ri)). The superscript pd stands for 
pair density. Similarly, the exact representation of E£ by 
integration over the non-linear adiabatic connection 

% = lj°° d tjj "(ri)ng(ri,r 2 ) ^^ dndr a , 

(15) 

where n|(ri,r2) is the correlation hole at interaction pa- 
rameter £, suggests a natural definition for the short- 
range local correlation energy per particle (see also 
Ref . El 




Note that even for the KS system (/i = 0), the local 
correlation energy per particle e^ =0:Pd (r) can depend on 
the adiabatic connection path followed. 

We have calculated these local energies per particle 
from accurate exchange and correlation holes associated 
to the system of Eq. using the accurate potentials of 
Sec. Inland the programs Molpro and CASDI |27j . 



and will not be shown here. We just mention that a 
careful comparison shows that the potentials with the er- 
fgau interaction are more short-ranged, confirming than 
the erfgau interaction achieves a better long-range /short- 
range separation than the erf interaction [TJ. 



III. SHORT-RANGE LOCAL 
EXCHANGE-CORRELATION ENERGIES PER 
PARTICLE 



The short-range local exchange-correlation energy per 
particle e^ c (r) associated to the exact E£ c is not uniquely 
defined. However, one can define and calculate some 
"physically realistic" £^ c (r) which can be compared to 
approximations like the LDA. We discuss now two such 
possible choices: local energies per particle obtained from 
integration of pair densities along adiabatic connections 
and local energies per particle obtained directly from po- 
tentials. 



B. Short-range local exchange-correlation energies 
per particle from potentials 

We now define another short-range local exchange- 
correlation energy per particle e£' local (n) as a function 
of n, by requiring that it yields of course the exact short- 
range exchange-correlation energy 



E 11 = 



n(r)^ local (n(r))dr, 



(17) 



but also the exact short-range exchange-correlation po- 
tential 



'd (ne^\n)Y 



dn 



5 (r) + C. (18) 



n— n(r) 



The potential is defined only up to an additive constant 
C and v£ c (r) is the exchange-correlation potential which 
goes to at infinity (calculated in Sec. ITTjl . 

For a given system, e"^ local (n) defines an exact local 
short-range exchange-correlation functional in the sense 
that it gives the exact exchange-correlation potential and 
energy. This approach has already been applied for 
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the Kohn-Sham scheme |2g; we briefly recall here how 
£xc° cai ( n ) is calculated. 

The condition of Eq. (|18|l implies 

V (ns^\n)) = ggg- WJ Vn = « c (r)+C)Vn, 



(19) 

« c M + c). (20) 



rf/i 

which becomes for spherically symmetric systems 

d(n(r)e^\r)) _ dn{r) 
dr dr 

Integration of Eq. I|20|) leads to 

n{r)e^\r) - (n(r)^ local (r)) ^ = 

^K c (r') + C)dr'. (21) 

To avoid divergence of the integral in Eq. (fTTj) . we must 
have (n(r)£^ local (r)) r _ +oo = 0; e^ local (0 is then ex- 
pressed as 



e^\r) = 



-1 



'' n ^ c (r')dr' + C, (22) 



n(r) J r dr' 
where the constant C is fixed by requiring that 



n{r)e^ oca \r)Airr 2 dr 

An f~> d(n(r)e^ oca \r)) ^ 
3 Jo dr 

An rdn{r) <Ar y dr + CN ^ ^ 



dr 



3 J dr 

with the number of electrons N = J n(r)Anr 2 dr . Sim- 
ilarly, using only the exchange or correlation parts of 
v% c (r) and Eg c in Eqs. and leads to the local 
exchange energy per particle e^ ,local (r) and local correla- 
tion energy per particle £^ ,local (r), respectively. 

We have calculated these local energies per particle us- 
ing the accurate potentials of Sec. Inland accurate values 
of Eg and Eg obtained from of Eq. fTTft . 



C. Results 

For monotonically decaying spherical densities, the 
map r — > n(r) can be inverted and the local exchange- 
correlation energy per particle can therefore be expressed 
as a function of the density n or, equivalently, as a func- 
tion of l/r s = (47T7J/3) 1 / 3 . The reason to use l/r s is 
that the exchange energy per particle of the Coulombic 
uniform electron gas g^ =0 ' unif is proportional to l/r s . 

For the He atom with the erf interaction, the two ac- 
curate short-range local exchange energies per particle 
eg' pd (l/r s ) and e^ local (l/r s ) are compared in figure |U 
together with the LDA local energy per particle. The dif- 
ference between the two accurate local energies per par- 
ticle, both giving the same correct short-range exchange 



energy Eg, illustrates the arbitrariness in the definition 
of local quantities from global ones. However, the differ- 
ences soften when fi is increased, i.e. when long-range in- 
teractions are removed. The LDA local exchange energy 
per particle is significantly different from both accurate 
local exchange energies per particle at the KS end of the 
adiabatic connection (/i = 0), but as (x increases the LDA 
local energy per particle better and better agrees with 
the accurate local energies per particle. At low densities 
(l/r s < 0.5), the LDA local energy per particle is close 
to the accurate eg' pd (l/r s ) while at high densities (for 
l/r s > 1.5) the LDA local energy per particle is closer to 
the accurate e£< local (l/r s ). 

In figure again for the He atom with the erf interac- 
tion, the accurate local correlation energies per particle 
e%' pd (l/r s ) and e^' local (l/r s ) are represented and com- 
pared to the LDA. The two accurate local energies per 
particle give the same correct short-range correlation en- 
ergy Eg but have rather different shapes, although the 
differences partly soften when /i increases. The LDA gen- 
erally overestimates the accurate local correlation ener- 
gies per particle at /j = 0, but again when /i increases 
the LDA is improved on average. The LDA local energy 
per particle is much more comparable to e^ ,pd (l/r s ) than 
e£' local (l/r s ) especially for large fj, and at low densities. 

In general, similar behaviors are obtained with the erf- 
gau interaction. However, the attractive character of this 
interaction for small values of \x (see Ref. ^5) can lead to 
some significant differences compared to the erf inter- 
action. Figure shows these differences in the accurate 
local correlation energy per particle £^ ,pd (l/r s ) at fi = 0. 
In comparison to the erf interaction, the local energy per 
particle obtained with the erfgau interaction exhibits in 
particular a positive contribution at very low densities. 
This emphasizes now the arbitrariness of the definition 
eg' pd (l/r s ) with respect to the choice of the adiabatic 
connection along which the pair density is integrated. 

Figure [7\ show the evolution of the local exchange- 
correlation energies per particle with the erf interaction 
for the Be atom. In comparison to the He atom, both 
e^ pd (l/r s ) and e^ local (l/r s ) exhibit a shell structure for 
small values of fj,, the contributions from the valence re- 
gion (r s < 0.5) and the core region (r s > 0.5) being easily 
identifiable. The LDA local exchange-correlation energy 
per particle differs significantly from the two accurate 
ones, especially in the core region which constitutes the 
major contribution to the exchange-correlation energy. 
Also, in this representation in term of l/r s , the LDA 
does not reproduce at all the intershcll jump. However, 
when /i increases, the contribution from the valence re- 
gion is progressively cut off and the shell pattern finally 
disappears. For large fj,, when only the contribution from 
the core remains, the LDA local energy is close to the ac- 
curate local energies. 

Finally, the accurate local correlation energy per par- 
ticle e^ ,pd (l/r s ) for the Be atom is reported in figure 00 
together with the LDA. The shell structure in this plot is 
even clearer. At n = 0, the LDA overestimates the accu- 
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FIG. 4: Short-range local exchange energies per particle with 
respect to l/r s for He with the erf interaction for fj, = 0, 1, 3: 
-£,lda (Jong.daghed curve) is compared to the accurate £j' pd 
[Eq. JT1J, solid curve] and e£' local [Sec. flTTBl . short-dashed 
curve] . 



rate local correlation energy per particle at all densities. 
When /i increases, the LDA starts to better reproduce 
the valence region but still overestimate the accurate lo- 
cal energy in the core region. Keeping on increasing fi, 
the contribution from the valence region vanishes, and 
the LDA starts to reproduce well the core region. 



IV. CONCLUDING REMARKS 

Non-linear adiabatic connections have been used to de- 
compose the energy of an electronic system into a long- 
range wave function part and short-range density func- 
tional part, the position of the frontier of this decompo- 
sition being determined by a single parameter \x. 

When the interaction range of the short-range energy 
component is reduced, corresponding to a rise of /i, the 
short-range LDA becomes an increasingly accurate ap- 
proximation to the short-range exchange-correlation den- 
sity functional. Indeed, as it has been verified in this work 
for a few atomic systems, the short-range LDA gives, 
for a sufficiently large good short-range exchange- 
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FIG. 5: Short-range local correlation energies per particle 
with respect to l/r s for He with the erf interaction for /i = 0, 
1, 3. £c' LDA (long-dashed curve) is compared to the accurate 
£^ pd [Eq. GS), solid curve] and e£,iocal [ Sec fTTTrl . short- 
dashed curve]. 
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FIG. 6: Accurate local short-range correlation energies per 
particle e^' pd [Eq. Iltjl l with respect to l/r a for He with the 
erf (solid curve) and erfgau (dotted curve) interactions for 
fi = 0. 



correlation potentials and local exchange-correlation en- 
ergies per particle. 

More formally, one can define an optimal local interac- 
tion parameter, /i opt (r), so that for /i > /x opt (r) the LDA 
reproduces, at a given precision, the exact short-range 
potentials or local energies at the considered point of 
space. Physically, l/^ op t(r) represents the maximum lo- 
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FIG. 7: Short-range local exchange-correlation energies per 
particle with respect to l/r s for Be with the erf interaction for 
= 0, 1, 3: e£c LDA [Eq. 10, long-dashed curve] is compared 
to the accurate e^ pd [Sec. 1111 All , solid curve] and e^ local 
[Eq. short-dashed curve]. 



FIG. 8: Short-range local correlation energies per particle 
with respect to l/r s for Be with the erf interaction for /i = 0, 
1, 3: £c' LDA (long-dashed curve) is compared to the accurate 
e£' pd [Eq. GHJ, solid curve]. 



cal interaction range over which the exchange-correlation 
effects are well transferable from the uniform electron 
gas to the inhomogeneous system of interest. In general, 
the larger the density at point r is, the larger /i op t(r) is. 
However, it is not conceivable in the method to use a in- 
dependent interaction parameter /i because of the wave 
function part of the calculation. 

In a simple system like the He atom, the improvement 
of the LDA upon increasing fi is quite uniform in space 
and one can basically define an global average optimal 
interaction parameter, flopt, and use it to fix the frontier 
of the decomposition, i.e. \i = ftopt- 

More generally, in strongest inhomogeneous systems 
like the Be atom, at least two regions of space are iden- 
tifiable, the core and valence shells, and it is meaningful 
to define separate averages of // pt(r) over each regions, 
/x™t e and rjop t cncc . As the density is higher in the core 
than in the valence, we have £topt e > A*opt ence - 

If one fixes the frontier of the long-range/short-range 
decomposition at the optimal value for the core shell, i.e. 
/i = /iopt e , the short-range functional part of the calcula- 
tion can be well treated in LDA but of course the comple- 



ment part of the core region and the entire valence region 
are assigned to the wave function part which can result 
in the need of a long expansion into Slater determinants. 

On the contrary, if one fixes the frontier of the decom- 
position at the optimal value for the valence shell, i.e. 
M = Mopt cnc °j the short-range LDA functional is not en- 
able to treat well the core region but can describe well 
that part of the valence region assigned to it. Moreover, 
the wave function part of the calculation has only to han- 
dle with the remaining part of the valence region which 
allows in general the use a short expansion into Slater 
determinants. 

Obviously, one would like to re-conciliate these two al- 
ternatives. For practical applications of the method to 
molecular systems of chemical interests, it is desirable 
to keep the wave function expansion minimal and thus 
to chose the frontier of the decomposition inside the va- 
lence region, i.e. /i = /i™ 1 ™ 00 . For properties essentially 
depending on the valence shell only, such as atomization 
energies, the error of the LDA in the core shell can then 
be ignored, or the core electrons can be replaced by a 
pseudo-potential. Otherwise, one must then find better 
approximations to the short-range exchange-correlation 
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functional which extend the domain of accuracy of the 
LDA toward larger interaction ranges (small fi). We hope 
that the presented results will help the construction of 
such approximations. 

Data for the atomic systems presented in this paper 
(and other systems) are available from the authors upon 
request. 
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